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ABSTRACT 

With an extinction-corrected V-band peak absolute magnitude of about —22, supernova (SN) 2006gy is 
probably the brightest SN ever observed. We report on multi-wavelength observations of this SN and its 
environment. Our spectroscopy shows an Ha emission line as well as absorption features which may be 
identified as Si II lines at low expansion velocity. The high peak luminosity, the slow rise to maximum, and 
the narrow Ha line are similar to those observed in hybrid type-Ia/IIn (also called Ha) SNe. The host galaxy, 
NGC 1260, is dominated by an old stellar population with solar metallicity. However, our high resolution 
adaptive optics images reveal a dust lane in this galaxy, and there appears to be an H II region in the vicinity of 
the SN. The extra-ordinarily large peak luminosity, ~3 x 10 44 erg s , demands a dense circum-stellar medium, 
regardless of the mass of the progenitor star. The inferred mass loss rate of the progenitor is ~ 0. 1 M yr~ 1 over 
a period of ~ lOyr prior to explosion. Such an high mass-loss rate may be the result of a binary star common 
envelope ejection. The total radiated energy in the first two months is about 1.1 x 10 erg, which is only a 
factor of two less than that available from a super-Chandrasekhar la explosion. Therefore, given the presence 
of a star forming region in the vicinity of the SN and the high energy requirements, a plausible scenario is that 
SN 2006gy is related to the death of a massive star (e.g., pair production SN). 

Subject headings: supernovae: general - supernovae: individual (SN 2006gy) - galaxies: individual 
(NGC 1260) 



1. INTRODUCTION 

SN2006gy was discovered by the ROTSE-IIIb telescope 
at the McDonald Observatory on UT 2006 September 18.3 
(Quimby 2006). The supernova (SN) was initially reported 
2" off the center of NGC 1260. Harutyunyan et al. (2006) 
obtained a spectrum on UT 2006 September 26 and reported 
a three-component Ha emission line: an unresolved nar- 
row line; an intermediate component with Full Width at Half 
Maximum (FWHM) of 2500 kms -1 ; and a component with 
FWHM of 9500 kms -1 . They suggested that the event was a 
type II SN. 

Prieto et al. (2006) reported that a spectrum of the SN, ob- 
tained eight days after the discovery, was suggestive of a dust- 
extinguished type-IIn event. However, the Balmer lines were 
symmetric, which is unusual for SNe in their early phases. 
Moreover, after correcting for two magnitudes of extinction 
(based on observed Nal lines), the absolute magnitude is 
about —22. They further reported that the position of the 
SN is consistent with the center of the galaxy, suggesting it 
is more consistent with an eruption of the active galactic nu- 
cleus (AGN) of NGC 1260. Foley et al. (2006) noted that the 
SN is offset by about 1" from the nucleus of NGC 1260. This 
fact along with a spectrum obtained six days after discovery, 
led them to suggest that SN 2006gy was a type Iln event. 

Here we report on multi-wavelength observations of 
SN2006gy. An independent contemporary analysis is pre- 
sented by Smith et al. (2006). 
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FIG. 1 . — r-band (filled circles) and i'-band (empty squares) light curves of 
SN 2006gy. Errors include the uncertainty in the absolute calibration. The 
bottom horizontal axis shows the time relative to the discovery of SN 20()6gy 
(2006 September 18). 



We initiated a photometric (g,r,i,z) monitoring program 
with the Palomar 60-inch robotic telescope (Cenko et al. 
2006). Not possessing pre-explosion images of NGC 1260, 
which are essential for accurate subtraction of the light from 
the host galaxy, we used archival R- and /-band images ob- 
tained with the Jacobus Kapteyn Telescope 5 on 1996 Jan- 
uary 13 and 1991 December 1, respectively. The r- and i'- 
band measurements and errors, presented in Fig. [T] and Ta- 
ble[TJwere produced by image subtraction using the Common 
Point-spread-function Method (CPM; Gal- Yam et al. 2004). 

On UT 2006 September 26 and December 18 and 19 we 
obtained spectra using the Low Resolution Imaging Spectro- 
graph mounted on the Keck-I 10-m telescope (LRIS; Oke et 

5 ING archive: http://casu.ast.cam.ac.uk/casuadc/archives/ingarch 
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TABLE 1 

Log of observations and measurements 



Date 


Tel 


: a — 

Magnitude 


Date 


; 

Tel 


: — 

Magnitude 


U r 2006 


& Band 




UI 2006 


& Band 


or Flux" 


09-24.9 


P60 b r 


14.80 ±0.10 


10-30.7 


P60 i 


14.29 ±0.12 


09-25.9 




14.78 ±0.10 


10-31.7 




14.30±0.12 


09-26.7 




14.73 ±0.11 


11-01.7 




14.31±0.12 


10-05.9 




14.45±0.10 


11-02.8 




14.22±0.12 


10-29.7 




14.32±0.10 


11-08.8 




14.30±0.12 


10-30.7 




14.28±0.10 


11-09.8 




14.34±0.12 


11-09.8 




14.40±0.10 


11-12.0 




14.30±0.12 


11-12.0 




14.33±0.11 


11-01.3 


P200 c J 


12.96±0.14 


09-24.9 


P60; 


14.90±0.12 


11-01.3 


P200 K s 


12.59±0.17 


09-25.9 




14.86±0.12 








09-26.7 




14.69 ±0.17 








10-05.9 




14.53±0.12 


11-20.4 


VLAX d 


186±80^Jy 


10-26.7 




14.30 ±0.12 


11-20.4 


VLA K d 


59±110/jjy 


10-29.7 




14.26 ±0.12 


11-23.2 


VLA Q d 


56±120AtJy 



a Observed magnitude or flux density of the SN. Magnitude errors include the 
uncertainty in absolute calibration, which dominates the errors. To convert 
specific-flux errors to 3-0" upper limits multiply the errors by 3. b Palomar 
60-inch (P60) magnitudes are given in the AB magnitude system. Abso- 
lute calibration was performed by fitting the Hipparcos BjVj and 2MASS 
(Skrutskie et al. 2006) JHK magnitudes of three nearby Tycho-2 (H0g et 
al. 2000) reference stars to synthetic photometry of stellar spectral tem- 
plates (Pickles 1998) in the same bands. The best fit spectral template 
of each star was used to calculate its synthetic magnitudes in the r- and 
('-bands. The uncertainty in this calibration process, calculated from the 
scatter between the zero-points derived from each Tycho-2 star, is about 
0. lmag. c Palomar 200-inch (P200) IR Vega-based PSF-fitting magnitudes, 
relative to the 2MASS star 03172629+4124103 within the field, as mea- 
sured with IRAF/DAOPHOT. d Center frequency of VLA bands are as fol- 
lows: 8.4 GHz (X), 22.5 GHz (K) and 43.3 GHz (Q). 

al. 1995). Spectra were also obtained on UT 2006 October 
28.3, 29.4 and November 25.2, using the Double Beam Spec- 
trograph (DBSP) mounted on the Hale 5-m telescope. The 
spectra are displayed in Fig. [2] 

On UT 2006 Nov 1 .3 we observed the event with the Adap- 
tive Optics system (Troy et al. 2000) equipped with the Palo- 
mar High Angular Resolution Observer (Hay ward et al. 2001) 
camera mounted on the Hale 5-m telescope. We used the 
wavefront reconstruction algorithm - denominator-free cen- 
troiding and Bayesian reconstruction (Shelton 1997), which 
delivered X" s -band images with 0"1 FWHM and a Strehl 
ratio of ~ 15%. We obtained 660s and 300s images in 
the K s and /-bands, respectively, using the high-resolution 
mode (25maspix~') and a 240 s .Ks-band image using the 
low -resolution camera (40maspix _1 ). Each frame was flat- 
fielded, background subtracted, and repaired for bad pixels 
using custom PyRAF software 6 . 

The field of SN 2006gy/NGC 1260 was observed by the 
Swift X-Ray Telescope (XRT) on 2006 October 30 and the 
Chandra X-ray Observatory on 2004 December 23 and 2006 
November 14 7 . For the Swift observations, assuming a Galac- 
tic neutral Hydrogen column density Nh = 1.3 x 10 21 cm -2 
(Dickey & Lockman 1990), and a power-law spectrum with 
index 1.8, we set a 3- (7 upper limit for the flux in the 
0.2-10keVbandof < 1.8 x lO^ergs- 1 cm" 2 . The Chan- 
dra observations reveal a variable source at the position of 
the nucleus of NGC 1260. The spatial coincidence lead us to 
attribute this source to an active galactic nucleus. In order 

6 PyRAF is a product of Space Telescope Science Institute, which is oper- 
ated by AURA for NASA. 

7 This latest observation was conducted under Director's discretionary 
time (PI: Pooley). 
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FIG. 2. — The spectrum of SN2006gy after subtracting the scaled SO tem- 
plate, and correcting for Milky-Way and NGC 1260 extinction (black line; 
see text). The blue line shows the same spectrum after subtraction of a 
third-degree polynomial fitted to the spectrum. The magenta line shows the 
spectrum of the luminous type-la SN 1991T at nine days post peak bright- 
ness, after the same processing. The spectrum of SN 1991T was redshifted 
by ~ 8500km s -1 in order that the possible Si II features in both spectra 
coincide. A zoom-in on the Ha emission line as observed by LRIS and 
DBSP is shown in the upper inset. The lower inset shows a section of the 
2-dimensional Keck spectrum of SN 2006gy obtained under good seeing con- 
ditions on 2006 December 18. The extension in the spatial (vertical) direction 
is an Ha emission near the SN location (as previously reported by Smith et al. 
2006). Technical details: For the LRIS spectra, with integration time rang- 
ing from 600 to 2400 s, we employed the 1"5 slit with the 400/8500 grating 
blazed at 7550 A, and the 600/4000 grism, in the red and blue sides, respec- 
tively. One exception is the December 19 spectrum that was obtained using 
the high-resolution R1200/7500 grating centered on the Nal lines. The DBSP 
observations were obtained with a 1"5 slit and R158/7500 and B600/4000 
gratings on the red and blue arms, respectively. The spectrum marked with 
October 28.9 is the sum of four spectra (total integration time of 1500s) ob- 
tained during the October ran. The integration time for the November DBSP 
spectrum was 900 s. 
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FIG. 3. — Left: K s band image. The SN (marked) is clearly resolved 
from the galaxy nucleus, right: /-band image after subtracting the best fit 
Sersic profile from the galaxy and a Gaussian profile from the SN using Gal- 
Fit (Peng et al. 2002). The Sersic model parameters are as follows: in- 
dex of 3.7, an effective radius of 34", an axial ratio of 0.51, and a position 
angle of 80deg. A dust lane (white band) is seen southward of the galaxy 
nucleus. Based on three 2MASS sources, we derived the galaxy nucleus po- 
sition: a = 03 ft 17 m 27/241, S = +41 24 m 18!'55 and the SN position (end 
numbers): 27*. 158 (a) and 18"88 (<5). The astrometric solution has rms of 
0"04 and 0"01 in a and S, respectively. 

to constrain the X-ray luminosity of the SN, we fitted the X- 
ray image with a model containing three components: A nar- 
row Gaussian centered on the galaxy position; a wide Gaus- 
sian centered on the galaxy position (i.e., diffuse emission); 
and a narrow Gaussian centered on the SN position. We find 
that the SN flux is consistent with zero, and that its flux is 
< 1.6 x 10 40 ergs _1 at the 3-C confidence level, assuming a 



SN 2006gy 



3 



photon index of 1.8, a distance of 73 Mpc to NGC 1260 and 
a neutral Hydrogen column density of Nh = 6.3 x 10 21 cm~ 2 
(in the Galaxy and NGC 1260). This field was also observed 
using the Swift Ultraviolet/Optical Telescope. We are await- 
ing late epoch observations in order to properly remove the 
host contamination. 

We performed radio observations of NGC 1260 with the 
Very Large Array (VLA) 8 on 2006 Nov 20 and 23 UT. The 
observations were obtained in continuum mode with a band- 
width of 2 x 50 MHz. We observed 3C48 (J0137+331) 
for flux calibration, while phase referencing was performed 
against J03 19+415. The data were reduced using the Astro- 
nomical Image Processing System. We did not detect a source 
at the position of the SN (see Table[TJ. 

3. SPECTRAL ANALYSIS 

In the optical spectra we identify two Na I absorption 
lines, one of Galactic origin and the other at the redshift 
of NGC 1260 (z = 0.019). Based on the ratio between the 
equivalent widths of these two absorption lines (e.g., Mu- 
nari & Zwitter 1997), we estimate the total extinction to- 
ward SN2006gy to be ~ 4.4 times the Galactic extinction 
(E B - V = 0.16; Schlegel et al. 1998), which gives E B -v « 0.7. 
The high-resolution spectrum, obtained on 2006 December 
19, resolves the Na I doublet. Based on this spectrum, we 
find that both the Galactic and NGC 1260 doublets have sim- 
ilar line ratios and are not saturated. We note that using the 
Na I-extinction correlation derived by Turatto et al. (2003), 
we find a total Eg-v extinction in the range 1 to 3.5 mag. 
Moreover, the extinction is derived by assuming that all the 
Na I absorption is of light emitted by the SN, rather than of 
light emitted by the host galaxy. Therefore, the extinction 
toward SN 2006gy is uncertain and this issue will require fur- 
ther study. 

In Fig. |2] we show the DBSP spectrum of SN 2006gy (black 
line) after the subtraction of a scaled SO galaxy template (Kin- 
ney et al. 1996). The template was reddened to account for 
Galactic extinction in the direction of the SN, and scaled so 
that the synthetic r — i color of the host-subtracted SN spec- 
trum matches the photometrically observed value at the same 
epoch. Next, we flux-calibrated the spectrum by requiring its 
r-band synthetic photometry to equal the observed magnitude 
of the SN at the same epoch. Finally, we corrected the spec- 
trum for extinction, assuming Eg-v ~ 0.7 mag and Ry = 3.08 
(Cardelli, Clayton, & Mathis 1989). From the final spectrum 
we find, at maximum light, an extinction corrected synthetic 
V-band magnitude (Vega system) of about 12.4 mag. 

Our spectra show an Ha and Hj3 emission lines with a P- 
cygni profile, characteristic of type-IIn SNe. We note that the 
equivalent width of the Ha line is decreasing with time. Inter- 
estingly, we detected several absorption features which may 
be Si n, S II, Fe II, Fe III and Ca II lines. Such lines are usu- 
ally observed in type-la SNe. However, we stress that the rela- 
tive lines strenghts and apparent low expansion velocities (i.e., 
~ 1000-2000kms _1 ) are peculiar, making line identifications 
tentative only. To emphasize this, we show in Fig.|2]the spec- 
trum of SN 1991T (Filippenko et al. 1992) at nine days from 
maximum light redshifted by 8500km s , and the spectrum 
of SN 2006gy at 42 days since discovery, after the subtraction 
of third-degree polynomials fitted to each spectrum. The lines 

8 The Very Large Array is operated by the National Radio Astronomy 
Observatory, a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



of SN 2006gy are narrower and red-shifted relative to these of 
SN 199 IT, indicating that SN2006gy had a lower expansion 
velocity compared to type-la SN. 

4. ENVIRONMENT 

NGC 1260 is an early-type galaxy within the Perseus cluster 
of galaxies. Its Heliocentric recession velocity is 5760 kms -1 
and its velocity dispersion is 201 ± 12kms _1 (Wegner et al. 
2003). Based on the recession velocity of the cluster the dis- 
tance modulus to NGC 1260 is 34.5 mag 9 . 

Our adaptive optics images (Fig.0 show that SN 2006gy is 
located 0"99 (projected distance 380pc), at a position angle of 
290 deg, from the nucleus of NGC 1260. A dust lane, passing 
about 300 pc (projected) from the SN location, is clearly seen 
in our galaxy-subtracted /-band image. Moreover, we confirm 
the detection by Smith et al. (2006) of an H II region in the 
SN vicinity (Fig.EJi. 

The Mg2 index of this galaxy was measured to be in 
the range 0.24-0.27 mag (Davis et al. 1987; Wegner et al. 
2003). This value, along with the synthetic spectral models of 
Vazdekis (1999), suggests that the metallicity of NGC 1260 is 
not low, [Fe/H]>-0.2. 

5. DISCUSSION 

With estimated peak absolute magnitude of V « —22, 
SN 2006gy is probably the brightest SN ever observed. The 
slow brightening, the peak luminosity and the Ha emis- 
sion line and the possible SN-Ia-like features suggest that 
SN 2006gy maybe related to the hybrid Iln/Ia SNe class (also 
known as type-IIa; Deng et al. 2004). The other possible 
members in the type-IIa group are SN 2002ic (Hamuy et al. 
2003) SN2005gj (Aldering et al. 2006); SN 1997cy (Ger- 
many et al. 2000) and SN 1999E (Rigon et al. 2003). 

Any model of SN 2006gy has to explain the spectral lines, 
the extra-ordinary peak luminosity of L p ~ 3 x 10 44 ergs _1 
(after correction for extinction), and a radiated energy over the 
first two months of E rai j ~ 1.1 x 10 51 erg (assuming 11,000-K 
black body which roughly matches the Rayleigh- Jeans slope 
in DBSP spectra). We note that even if the extinction in 
NGC 1260 was overestimated and the SN light suffers only 
from Galactic extinction, the total radiated energy within the 
first two months is about 3 x 10 50 erg. 

The high peak luminosity suggests that the blast wave from 
the explosion efficiently converts the mechanical energy to ra- 
diation. This mean the shock has to be radiative which re- 
quires the circum-stellar medium (CSM) density to exceed 
10 6 cm~ 3 . Moreover, the conversion of mechanical energy 
of an explosion to radiation requires that the ejecta sweep up 
matter with comparable mass. The slow rise time, t p ~ 50 d 
to peak luminosity implies that the dense region has a size of 
at least R ~ v s t p ~2x 10 15 cm, where v v ~5x 10 8 cms~' is 
the speed of the blast wave. The peak luminosity, L p , requires 
density of the order, n ~ L p / (2nR 2 v^) ~ 10 10 cm~ 3 . Assum- 
ing an upper limit on v s of 10 9 cms _1 at early times, the min- 
imum mass contained within this radius is J> O.2M . The 
gradual decrease in the radiated energy, and possibly lower 
expansion velocities would easily bring it closer to a solar 
mass. The mass loss rate by the progenitor has to be stu- 
pendous, M ~ lM Q /(t p v s /v w ) ~ 10 _1 M Q yr _1 , over a time 
scale of at least about lOyr, where v w — 200kms~ 1 is the 

9 Assuming Hubble parameter Hq = 72 km s Mpc~' , matter content 
Q m = 0.27, and dark-energy content of Cl\ = 0.73. 
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speed of the progenitor wind (Smith et al. 2006; v w ~ 130— 
260km s -1 ). Finally, the high CSM density accounts for the 
lack of substantial X-ray and radio emission (being absorbed 
by photoelectric and free-free absorption, respectively). 

SN2006gy shares some properties with type-IIa and type 
Iln SNe. Type-IIn SNe are most plausibly the result of a core 
collapse SN embedded in dense CSM, while Ha events have 
been explained as thermo-nuclear explosions taking place in 
a dense medium (e.g., Livio & Riess 2003; Han & Podsiad- 
lowski 2006). The thermo-nuclear model is attractive from 
a spectroscopic perspective. In the context of type-la SNe, a 
possible explanation to the high-mass loss rate is that it is the 
result of a common-envelope phase in a binary system (e.g., 
Taam & Ricker 2006 and references therein). This scenario 
was suggested by Livio & Riess (2003) to explain the prop- 
erties of SN 2002ic, and is consistent with the inferred high 
mass loss rate and its velocity (i.e., ~ 200kms~'). However, 
this scenario requires the ejection of matter from the progeni- 
tor to shortly precede the SN explosion (Chugai & Yungelson 
2004). Moreover, the total kinetic energy of la events is lim- 
ited to about 1-2 x 10 51 erg (Khokhlov et al. 1993), and it 
can get up to 2.5 x 10 51 erg for super-Chandrasekhar models 
(cf. Yoon & Langer 2004). Therefore, unless we consider- 
ably over-estimated the extinction, the total radiated energy 
of SN 2006gy in the first two months alone is challenging for 
type-la-like SN models. 

Smith et al. (2006), noting that the envelope of a massive 
star (> 100M©) contains a reservoir of thermal energy that 
can power the SN, suggested that such a star was the pro- 
genitor of SN 2006gy. However, most of the thermal energy 
will be lost due to expansion and the ability of the photons to 
leak out is limited by the long diffusion timescale for photons 
(;>monfhs; e.g., Kulkarni 2005). Therefore, it will be diffi- 
cult for this specific model (alone) to explain the high peak 



luminosity of SN 2006gy. 

Along the general lines of previous suggestion by Benetti 
et al. (2006; for SN2002ic) and Smith et al. (2006; 
for SN2006gy), we speculate that the large energy budget 
for SN2006gy may hint at a highly energetic explosion (~ 
10 52 erg), from a massive stellar progenitor. Two possibili- 
ties are an CSM-embedded collapsar (e.g., Woosley & Mac- 
Fadyen 1999), or a pair production SN (e.g., Ober et al. 1983; 
Smith et al. 2006). Pair production SNe, however, require 
low-metallicity progenitors, but it may be possible to over- 
come this requirement by the merger of two massive stars. 
We further speculate that such a merger may be responsible 
to the high mass-loss rate (e.g., common envelope ejection). 
SN 2006gy and other Ha (and maybe many Iln) events may 
result from one of these energetic explosions that are able to 
produce ~ 10 52 erg. 

The general issues of the large energy release into a dense 
CSM have been discussed for some type Iln events (e.g., 
Chugai et al. 2004, Gal-Yam et al. 2006). For reasons we 
do not understand the explosion is preceded by a phase of stu- 
pendous mass loss. The mass and geometry of the hydrogen 
envelope may determines the outcome of the explosion (e.g., 
Ha or Iln). Finally, the rarity of such energetic SNe reflect the 
rarity of the progenitors. 
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